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' We consider maximally supersymmetric U{N) Yang-Mills in (1 + p)-dimensions for p < 3. In 

Q the 't Hooft large N limit this is conjectured to be dual to N Dp-branes in the decoupling 

r^i limit. At low temperatures T <^ ^i/(3-p) governed by the dimensionful 't Hooft coupling 

, A, supergravity black holes predict the free energy density goes as ~ NT (^^p) and the 

J> expectation value of the scalars goes as ~ T^-p , with dimensions made up by A. The purpose 

^^ of this work is to explain the origin of these peculiar powers of temperature. We argue that 

0^ these powers naturally arise by requiring that the low energy moduli of the theory become 

cn 



strongly coupled at low temperature. As an application, we consider the BMN quantum 
mechanics that results from a supersymmetric deformation of the p = theory. The black 
C<^ holes dual to this deformed theory have not yet been constructed, and our analysis can be 

used to make an explicit prediction for their thermodynamic behaviour. 
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1. Introduction 

Maldacena's AdS-CFT duality [1] has remarkable potential to allow us to understand the 
most basic questions in quantum black hole physics. The simplest realisations from the field 
theory side of the correspondence are those arising from maximally supersymmetric Yang- 
Mills (MSYM) in (1 + p)-dimensions which is dual to the decoupling limit of Dp-branes [2], 
the original AdS-CFT case being p = 3. Whilst these field theories are straightforward to 
write down and are well defined for p < 3, it has so far not been possible to solve them. In 
the cases p < 3 the gauge theory is strongly coupled in the IR and for the conformal p = 3 
case, the dimensionless coupling must be taken to be large in order to achieve a weakly curved 
gravity dual. The challenge is to directly solve the thermal physics of these strongly coupled 
gauge theories and extract the dual quantum behaviour of black holes. ^ 

A fascinating mean field approach was developed by Kabat, Lifschytz and Lowe [5-7] for 
the p = case to directly solve the thermal theory, and extract the thermodynamics, with 
some success. Very recently there has been exciting renewed efforts to develop this method 
further [8] . Beyond these works there has been a distinct lack of analytic progress on explicit 



^For D1-D5 brane systems one obtains a realisation of AdS-CFT where the CFT is (1 + l)-dimensionaI [3]. 
In these cases the entropy of the dual BTZ black holes can be exactly reproduced using CFT methods [4] . It 
is worth emphasising that whilst this is a great success, the CFT is not directly solved. It is only the entropy 
that can be computed, and other quantities cannot be solved for. 



calculation of the thermal theory for any p < 3. Over the past decade an exciting new 
direction has been the direct numerical Monte Carlo simulation of these gauge theories at 
finite temperature on the lattice (or 'non-lattice') in the large N limit relevant for Maldacena 
decoupling. While the conformal p = 3 case is perhaps the most elegant to consider from an 
analytic perspective, from the gravitational point of view the quantum black hole physics is 
essentially the same in all the cases of p < 3. For the numerical approach lower dimensions is 
a big advantage, meaning less degrees of freedom to simulate and a simpler continuum limit. 
Indeed for p = 0, the BFSS quantum mechanics [9], the gauge theory is a quantum mechanics, 
hugely simplifying the continuum limit. In a lattice approach this implies that complicated 
supersymmetric lattice actions are not required [10] and an elegant and powerful non-lattice 
or spectral approach [11] can be taken since the gauge can be entirely fixed up to the Polyakov 
loop freedom. The first simulations of this p = have directly seen thermodynamic behaviour 
consistent with the dual predictions from supergravity [12-15], and have even allowed the a' 
corrections to supergravity to be extracted [16]. The BMN [17] deformation of the p = 
quantum mechanics, which we discuss later, has been studied [18]. Going beyond p = 0, 
the first thermal simulations of the p = 1 case on a spatial circle [19] have seen evidence 
of an expected Gregory-Laflamme phase transition [20-24]. The higher p cases are likely to 
be practical in the future, with the first work appearing recently on the conformal p = 3 
case [25]. For excellent recent reviews on this direct Monte Carlo approach see [26,27].^ 
Beyond equilibrium thermodynamics an exciting numerical approach has been employed [33- 
35] to study dynamics of initial states that thermalize and to understand how information is 
'scrambled' [36]. 

In the case of p = 3 the conformal invariance specifies that energy density e ~ N'^T^ for 
temperature T. For the cases of p < 3 there is more texture in the thermodynamics.^ Lacking 
conformal symmetry, the 't Hooft coupling A = Ngyj^ is dimensionful and one may form the 
dimensionless temperature t = T /\^'^^~p> . Now the energy density may have an interesting 
dependence on this dimensionless parameter. For p < 3 and high temperatures, t ^ 1, 
these theories are weakly coupled, and hence have energy density going as, e ~ N'^T^^^, 
being simply determined by the number of UV degrees of freedom in the theory. However, 
in the opposite low temperature limit, where as A^ — )■ oo then t is finite but small, so t <C 1, 
the theory is strongly coupled and a dual supergravity description in terms of a Dp-brane 
black hole becomes valid, predicting a particular N dependence and power law scaling with 



^ A number of other numerical works have simulated supersymmetric gauge theories to test holography and 
Matrix theory in the zero temperature context [28-32]. 

^In fact for p < 3 the canonical ensemble does not exist in these cases as the Dp-branes may Hawking 
evaporate from the decoupling region. This Hawking radiation rate is subdominant in A'^, and cannot be seen 
directly in the supergravity Dp-brane black hole thermodynamics, and at large N we may think of there being 
a metastable thermal state. More formally, addition of a finite but small mass term cures this IR divergence 
of the canonical partition function, and in the limit of small mass will not change the leading N behaviour. 
For an extended discussion of this in the context of p = see [15] where the consequences for Monte-Carlo 
simulation are discussed. For a recent estimate of this Hawking radiation rate in the string dual see [8] . 



temperature as, 

2(7-p) 1+p 

e~ Af^t (5-rt A3-P , t<l. (1.1) 

It is then a challenge for direct tests, numerical or otherwise, to reproduce both the correct 
temperature power law dependence, and also the overall coefficient. As mentioned above in 
the case of p = this has been achieved with great success using the direct numerical Monte 
Carlo approach. 

One might hope that just as for the conformal p = 3 case, the power law dependence of 
energy density with temperature can be analytically understood directly from gauge theory, 
even if calculating the precise coefficient cannot be achieved. It is the purpose of this paper 
to confirm this expectation. Our principle motivation is that so far the direct approaches 
(mean field or Monte Carlo) to solving p < 3 have not assumed this power law dependence, 
but rather have tried to reproduce it. However, clearly if it has a simple analytic origin in the 
gauge theory, one might hope to significantly improve these direct approaches by building in 
this analytic knowledge. For example, in the numerical approach one might hope this would 
yield understanding of the most relevant field configurations and could be used to improve 
the Monte Carlo sampling. For the mean field approach one might develop better informed 
Gaussian approximations to expand about. 

Our starting point is the work of Smilga, who in the case of p = considers the dynamics 
of the low energy moduli of the theory [37]. Smilga considers the quantum effective action 
of these moduli. The 1-loop correction to the moduli action is well known (the famous u ^/r^ 
term in Matrix theory [9,38]) and Smilga supposes that the configurations where the tree 
level and 1-loop terms are of the same order, ie. the moduli are strongly coupled, reproduce 
the correct temperature scaling. Making certain assumptions, the strongest being to assume 
the correct A^^ dependence, so taking e r^ N'^t'^X^'^ for some q, Smilga then deduces the 
correct power q = 14/5 that arises from the dual gravity by equating the leading tree level 
and 1-loop moduli action terms. Hence the peculiar dynamical scaling of energy density with 
temperature is given an origin within the MSYM. This approach is highly reminiscent of the 
work of Horowitz and Martinec [39], Li [40], and Banks, Fischler, Klebanov and Susskind [41] 
in the Matrix theory context to explain the parametric dependence of energy with horizon 
radius for Schwarzschild black holes using precisely the same logic, namely taking the moduli 
action to be strongly coupled by equating the tree and 1-loop terms in the moduli action. 

Motivated by these earlier Matrix theory works, and Smilga's application of these ideas 
to the Maldacena limit of DO-branes, we may ask various questions; 

• In the 't Hooft limit with finite but small t <^ 1 can we explain the temperature power 

2(7-p) 1+p 

law dependence e ~ N t (^-p) X^-p for all p < 3, not just p = 07 

• Can we explain it without assuming the correct N dependence (as Smilga assumed 
in [37])? 



• In addition to the energy density, Maldacena duality [2] predicts the behaviour of the 
scalar vevs to have a particular power law dependence on temperature. Can we predict 
this too? 

2(7-p) 1+p 

• The power law dependence e ~ N t (^-p) X'^^-v only arises when the dual supergravity 
is valid, at low temperatures t <^ \. Indeed, as we stated above, at high temperature 
e ~ A^^j^i+p fQ], ^]^g cases p < 3. Can we understand why this temperature scaling is 
only obtained in this particular low temperature limit (by which we mean t finite but 
small as A'' — )• cxd)? 

In this paper we shall address these questions. We will consider the thermal effective 
action for the moduli of the p < 3 theories. We will show that requiring the moduli action 
becomes strongly coupled, following the ideas of [37,39-41], naturally yields both the correct 
A^ and t dependence predicted by the dual Dp-brane black holes for both the energy density 
and the scalar expectation values for all p <2>. In addition we will explain why these t scalings 
only arise at low temperature where t is finite but t ^ 1. 

Beyond explaining the origin of the power law temperature scaling of the energy density 
and scalar vevs predicted by the usual dual Dp-brane black holes, these ideas may be used 
to make quantitative predictions for supergravity black hole solutions that have not yet been 
found. In particular in the case of p = MSYM one may deform this BFSS gauged quantum 
mechanics preserving the 16 supercharges and introduce a mass parameter /i, to give the 
BMN quantum mechanics [17]. As emphasised in [15] this is a particularly attractive system 
to study directly by Monte Carlo methods, as it is low dimensional, it has a IIA supergravity 
dual even with the mass deformation, and the addition of the mass naturally cures the IR 
singularities of the p < 3 theories hence rendering the canonical partition function finite. 
Despite this fact, there has been relatively little work on it, the first Monte Carlo study 
being [18]. The zero temperature IIA supergravity duals are known from the work of Lin [42] 
and Lin, Lunin and Maldacena [43], but finite temperature black holes with the correct UV 
asymptotics are not yet known, and this remains an important outstanding gravitational 
problem that must be solved to provide predictions for direct tests on the MSYM side. We 
argue that in the regime where supergravity is good, our simple analysis strongly constrains 
the dependence of the energy density and scalar vevs on the temperature and mass //. Since 
we do not yet know the gravity dual solutions, these constraints may be regarded as a gauge 
theory prediction for supergravity. 

The structure of this paper is as follows. In the next section we will briefiy review the 
supergravity predictions for the energy density and scalar behaviour at low temperature for 
(1 +p)-dimensional MSYM. Then in the following section 3 we will give the simple argument 
to explain the origin of the peculiar power law dependence on temperature in these quantities. 
We will consider the quantum corrections to the low energy moduli space of these theories, and 
requiring these to be strongly coupled at low temperature will reproduce the correct power 
laws. For readability, we will not discuss the somewhat technical derivation of these quantum 
corrections until the later section 5 of the paper. In section 4 we will discuss the BMN 



deformation of the p = theory, and using the same arguments, deduce constraints on the 
low temperature behaviour of the thermodynamics and scalar vevs for the, as yet unknown, 
supergravity black holes. In section 5 we will give the deferred details of the calculations of 
the thermal quantum 1-loop corrections to the classical moduli action. We conclude with a 
brief discussion. 

2. MSYM and supergravity black holes 



We begin with maximally supersymmetric U{N) Yang-Mills theory with coupling qym- In 

Jym- 



the large A^ 't Hooft limit the natural coupling becomes A = Ngy.,. The bosonic part of the 



action is given as. 



Sb = ^ j dtdxP Tr -^F^, - ^D^$^Z)^$^ + J [<^' , <^'' 



(2.1) 



and the fermionic action is, 



Sf = -—JdtdxPTr^i^^'D^ - i [y^',-]) ^ (2.2) 

where x^ with /x = 0, . . . ,p are the world volume coordinates with time x^ = t and spatial 
coordinates x* with i = 1, . . . ,p. Then $ with / = 1 + p, . . . , 9 are the spacetime scalars 
representing the transverse degrees of freedom of the branes. They are N x N Hermitian 
matrices transforming in the adjoint of the gauge group. The fermion ^ is a (l-|-9)-dimensional 
Majorana-Weyl spinor and also transforms in the adjoint, and the collection {7'^,7 } are 
the appropriate 50(1,9) gamma matrices. As usual, F^j^ = df^Ai, — d^Af^ + i[^^, A^], and 
Dfj, = dfj, — i[A^, •], with A^ the gauge field potential transforming in the adjoint. 

Maldacena duality [2] states that at finite temperature there is a dual closed IIA (p even) 
or IIB (p odd) string theory description of this gauge theory in terms of the decoupling limit 
of thermal Dp-branes. In the large A^ 't Hooft limit this string theory description may reduce 
to one in IIA or IIB supergravity, where the string frame metric takes the form. 



7~p 



ds^ = a' -^={-fdt^ + dx'dx') + y^ {U~^^ + C/ Vci!^2 



7-p 



nU) = l-{'^] ' , a, = 2^--.^r(l?) (2.3) 

and the dilation goes as, 

e*=i(2.)->,)^(3^)""^ (2,4) 

and in addition the (1 -|-p)-form potential carries the charge of the A^ Dp-branes. Here t and 
x* are the Dp-brane (1 -|- p)-dimensional world volume coordinates. The radial coordinate 



U is identified with an energy scale associated to the expectation value of the scalars in the 
MSYM. The horizon is at t/ = Uq, implying the thermal scale for the scalars is C/q so that, 



Uo^^^iTr^'^') (2.5) 

at finite temperature. 

In order for the supergravity solution to consistently describe the full closed string theory 
physics we require that both stringy and a' corrections are small. Stringy corrections are 
characterised by the string coupling, the dilaton. The a' corrections are characterised by 
the size of the curvature of the solution relative to the string scale set by a' . The radius of 

3-p 

curvature R goes as, a' /B? ~ {U /\^'^^~'P'^ ^ . We will consider the large N 't Hooft limit, 
where it is natural to take, 



and we see from (2.4) that in this large N limit the solution is well described by semiclas- 
sical string theory since the dilaton will everywhere be small. However, in order to ensure 
that supergravity gives a good description we shall in addition require that the curvature a' 
corrections are small, and so we must take, 

[/ , [/o < A^/(3-p) ^ (2.7) 

We emphasise that we require these quantities to be small, but finite, in the large iV — )• oo 
limit. 

The expectation value for the energy density e of MSYM is given by the energy density 
of the Dp-brane solution above extremality, and yields. 



;^(i+p)/(3-p) "P^' \xyi^~p)J ' ^ 9-p 

as deduced either from the Dp-brane solutions before taking the decoupling limit, or else 
from holographic renormalization that has been developed for these non-conformal Dp-brane 
geometries [44,45]. Defining the dimensionless temperature t = T/X^'^^~P' then semiclassical 
black hole thermodynamics yields, 

Cpt5-p , cp= ^ . (2.9) 



The condition Uq ^ A^' ''^~*''' above for supergravity to be valid translates into the condition, 
i <C 1 on the dimensionless temperature. Thus we see that gravity makes the following 
parametric predictions at low temperature; 

e ^ NH^-& X^^+py^^-P^ , wl(Tr$/$/) ^ [/^ ^ t^A^/^^-^^ (2.10) 

for t < 1. 



3. The basic argument 

We wish to consider putting our MSYM at finite temperature T = 1//3. It is convenient then 
to use the Euchdean time formulation, where Euchdean time r = it is periodic with period 
/3, and the action S = iS with S the Euchdean action, and A^ = iA^ for the gauge field. 
Hence now x^ = (r, x*). It is also convenient to introduce vector notation for the SO{9 — p) 
symmetry index, so now <1> — t- $. 

A set of classical vacua of this theory are gauge equivalent to the configurations, 

Kb = "a'^afe , ^ab = (faSab (3.1) 

for real constants aa and (pa, with fermion fields being trivial, and a,b = 1,. ..,N are the 
colour indices of the adjoint hermitian matrices. Such a configuration breaks the U{N) gauge 
symmetry to U{1) . The constants aa then represent constant gauge field configurations for 
this U{1) gauge symmetry. The (l)a give the transverse displacement of the A^ Dp-branes. 
We may promote the constants Oa, (j)a to slowly varying real scalar moduli fields, aa{T,x), 
(j)a{T,x), and then these have a classical moduli space Euclidean action, 

S''''^'-'^'^' = y y drdxP Yl {l^^^- ■ ^f^^- + i^M-^a ') • (3-2) 

where F^j^^ = d^a^a — d^a^a are the field strengths for the A^ U{1) gauge fields. Considering 
non-zero configurations for these moduli fields implies that the off diagonal components of the 
bosonic matrices A^, $ and fermionic spinor matrices ^ are massive degrees of freedom. The 
ab components of these matrices have mass squared ~ |^a — 0fep. Considering configurations 
of the moduli such that all these masses are large, one may integrate out the off-diagonal 
degrees of freedom to yield quantum corrections to this classical moduli space action. Whilst 
with maximal supersymmetry there can be no corrections to the potential and kinetic terms, 
in fact one does find higher derivative quantum corrections. 

We now apply a simple extension of the standard calculations familiar from consideration 
of Matrix theory [9,38,46] to deduce various 1-loop quantum corrections to this leading tree 
level moduli space action. We will study these quantum corrections up to 4 derivatives in the 
scalars, and ignore terms involving derivatives of the gauge field moduli. At 1-loop the terms 
come in two varieties; those that are independent of temperature (and hence are generated 
at zero temperature), and those that explicitly depend on temperature. We denote these 
the non-thermal and thermal 1-loop corrections respectively. We will now state the relevant 
results here, and defer their calculation until the later section 5. 

Firstly we consider the non-thermal correction. Up to 4 derivatives in the scalars, the 



non-thermal 1-loop correction has a single contribution of the form, 



S. 



E,l~loop 
non— thermal 



drdxP ^ 



r(^) / {^^,(^a-4>b]■^4^a 



a^Ua-0b -^M^a 



a<h (47r) 2 



|7-p 



^y.y^a-^h)-d^[4>a-<i^h)) 



\\2 



\(pa - (t>b 



7-p 



+ ... 



where supersymmetry ensures that the potential and 2 derivative kinetic terms are uncor- 
rected. We see that there is no explicit temperature dependence beyond the period of Eu- 
clidean time, and that the contribution, being negative, is attractive. This term is the gener- 
alisation of the famous v^/r'^ correction familiar from Matrix theory [9]. The dots . . . above 
are terms involving higher derivatives of the scalars, which on dimensional grounds take the 
schematic form, ~ (^{dcf))^ / 4?^^^ ((9(/>)^/(/)'') , for integer n > 1 (see for example [47] for the 
case p = 0). As we shall see later, these will be suppressed relative to the leading 4 derivative 
term. In addition, recall that we are only keeping track of terms involving derivatives of the 
scalars, and we similarly expect 4-derivative F^/(/>^~^ corrections for the gauge fields [48,49]. 
In addition to these temperature independent terms, finite temperature breaks the su- 
persymmetry and generates corrections for all even derivative orders that explicitly depend 
on temperature. However, an important feature of the supersymmetry is that due to the 
boson-fermion cancellations, these thermal corrections are exponentially suppressed in the 
limit. 



(3.3) 



»1 . 



(3.4) 



As we shall later see, it is precisely this limit that we will be interested in. We find the leading 
1-loop thermal correction in this limit to be, 



qE,1—Ioop 
thermal 



a<b 



UaU* + UbU^ 



IJl+P 



ma-M /^i^ _ ^^1 



p/2 



(3.5) 
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1 



X I — - — ^ — — -|- higher derivative terms | x 1 -|- O — ^ =^ 

(2vr)p/2 ; I \ma-M 



' = e* * " are the Polyakov loops around the Euclidean time circle. In 
b\ ^ I the two derivative corrections, being exponentially suppressed, are 



where Ua = e* 
this limit f3\(f)a 

subdominant to the tree level kinetic two derivative terms in (3.2), and the four derivative 
corrections are subdominant to the non-thermal 1-loop four derivative term of (3.3). As we 
shall therefore not require details of these terms, for simplicity we have suppressed them 
in the above expression, giving only the zero derivative potential term which is the leading 
zero derivative correction to the moduli action. This 1-loop zero derivative correction, the 
thermally generated potential, is the generalisation of that computed by Ambjorn, Makeenko 



and Semenoff in the case of p = [46].^ In the later section 5 some of the higher derivative 
terms are calculated and their explicit form is given. 

We now use an analogous approach to [37,39-41] in order to deduce the A^ and T scaling 
of the quantities of interest. One important difference to these works is that we will consider 
the theory at finite temperature using the canonical Euclidean approach rather than a Hamil- 
tonian Lorentzian point of view. This was first used in the Matrix theory context in [50] and 
appears to have the virtue that the N scaling need not be assumed, but emerges naturally. 

Having given the 1-loop corrected moduli action, we now show that the peculiar para- 
metric temperature dependence predicted by supergravity emerges precisely from requiring 
that the loop expansion becomes strongly coupled. Firstly we estimate all the (f)a to be of the 
same order which we write as (j). Furthermore we estimate that the differences between these 
moduli are also of the same scale, so that, 

(t>a ^ (f>a - (t>b ^ 4> ■ (3.6) 

At finite temperature Euclidean time is periodic, setting a classical length scale /?. We now 
make our key assumption, namely that it is this classical scale that determines the behaviour 
of the derivatives, so that we estimate, 

1 

r 

We emphasise that from the perspective of the leading classical moduli action, the quantum 
length scale A~^' (^"^^ naively plays no role, whereas the classical thermal scale /3 does. Hence 
we believe this is the most natural estimate to make, but we note that the remainder of this 
argument hinges on this point. We now proceed to estimate when the loop expansion becomes 
strongly coupled, and do this by equating the classical moduli action for the scalars with the 
1-loop correction, ignoring coefficients that are order one. We assume that the condition in 
equation (3.4) holds, so that, 

/3|4 - A\ ~ /30 » 1 (3.8) 

and then the 1-loop corrections are dominated by the non-thermal terms, with the thermal 
terms being exponentially suppressed. In addition it is worth emphasising that the 4 derivative 
1-loop non-thermal term given in equation (3.3) also dominates all the higher derivative 1-loop 
non-thermal terms, since as discussed below that equation, these are related by increasing 
powers going as (dcp)"^ / (j)'^ ~ l/{l3(p)'^ and hence are power law suppressed relative to the 4 
derivative term in this limit. A key point is that we will shortly check the self-consistency 
of this assumption finding it holds precisely in the low temperature regime t <C 1 dual to 
supergravity. 



5^<^a ~ df,(l)a - df,(j)b ~ ^0 • (3.7) 



''We also note that as discussed in [15] for p = 0, the 1-loop approximation becomes better in the regime 
/3|<j!>a — </>b| 3> 1 where for such a potential the moduli are essentially massless, and this leads to a divergence in 
the canonical partition function (at the subleading 0{N) rather than 0{N^)) as for p < 3 one should perform 
an integration over these moduli fields (due to Coleman-Mermin- Wagner) . 



The classical moduli term is estimated as, 

where we use f dr ~ /3, J2a ~ ^ ^^^ d^'Pa ~ 4>/ P- We similarly estimate the 1-loop 
contribution which is dominated by the 4 derivative non-thermal term to obtain, 



SE,.-ioo, ^ slt-}Z.^^, ~ \drdx^ Y: 



<b 

iV2 






dx^ TTTT^ (3-10) 

where we have taken, X]a<6 ~ -^^- Now equating the classical and 1-loop terms yields the 
estimate, 

2 
gE,classical _ gE,l-loop ^ ^ _ ('/JAS^) '"'' A^^ = t^A^^ (3.11) 

for the scalars. Having made the assumption (3.8) which suppresses the 1-loop thermal 
corrections and higher derivative non-thermal corrections, we must now check it is consistent. 
Using our estimate, 

3-p 

K/3(/>~ (/JA^^V""" =t"5^ (3.12) 

and hence we see for p < 3 then our assumption is consistent when t <^ 1. Hence the above 
analysis is consistent precisely in the low temperature regime where supergravity black holes 
give the predictions we wish to derive. At higher temperatures one would have to consider 
the contribution from the 1-loop thermal terms and higher derivative non-thermal terms. 
The thermal vev of the scalars is then estimated as, 



< ^Tr$^$^ > = J<J^Y1 1*^"!' > ~ <^ ~ t-AA3-P (3.13) 

where again ^^ ~ N, and we recall that in the regime t <^ 1 where our analysis is self- 



consistent then supergravity predicts \/< -^Tr$^<I>^ > ~ t^-pX'^-p which exactly agrees with 
our estimate. Encouraged by this we proceed to consider the free energy. Now in this low 
temperature regime we may use (3.9) to estimate the classical and 1-loop terms which we 
have equated to give. 



2 „ r / 1 X -Sup 



dx^ j^cp^ ~ A^ y dxPN' (/3A=^)"^-'' . (3.14) 
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The free energy, F, is given as I3F = {S^). Hence we may make estimates of the free energy 
density /, and energy density e (which for power law dependence on temperature has the 
same parametric dependence as /); 

P n / 1 \-§3; o 2(7-p) 1+p 

/3/~ A3^ivM/3A3^j "^ =^ e , f r^ nH^^=^ X^^ . (3.15) 

This is precisely the parametric dependence predicted by the supergravity in the regime t ^ 1 
where our estimates are self-consistent. 

Thus the condition that the quantum corrections to the classical moduli space action be- 
come strongly coupled at finite temperature precisely reproduces the supergravity predictions 
for the energy density and scalar vev in the regime t <C 1. In this temperature range we see 
that the 1-loop quantum corrections are dominated by the 4-derivative non-thermal terms 
(those with no explicit temperature dependence, that also are present at zero temperature), 
with the thermal contributions (those with explicit temperature dependence) being exponen- 
tially suppressed relative to these, and the higher derivative non-thermal corrections being 
power law suppressed. At higher temperatures, i ~ 1, the thermal 1-loop corrections and 
higher derivative non-thermal 1-loop terms are not suppressed relative to the leading non- 
thermal correction, and presumably the condition that the 1-loop corrections become strongly 
coupled is then much more complicated, and the simple power law temperature dependence 
of energy density and scalar vev will no longer apply. This is seen for p = Monte-Carlo 
simulations [12,13]. 

In the next section we use the simple estimates above to make a prediction for the 
thermodynamic behaviour of certain black holes in IIA supergravity, those dual to the BMN 
quantum mechanics. Then the remainder of the paper, section 5, is spent computing the 
1-loop quantum corrections given above in equations (3.3) and (3.5). 

4. A prediction for black holes dual to the BMN quantum mechanics 

Now we focus on the p = case, the BFSS quantum mechanics [9], and add to the action a 
term Sbmn that includes mass and Myers terms [51] that break the global SO(9) symmetry 
to 50(3) X 5*0(6) whilst preserving the 16 supercharges of the theory [17]. This theory is the 
BMN quantum mechanics. The bosonic part of this is. 



^bosonic ^ N ,^^^ 



^BMN 



2A 



1=1 7=4 I.J,K=1 



(4.1) 



After this mass deformation of the BFSS theory, in the 't Hooft limit we have been considering 
there remains a IIA string dual description. The zero temperature IIA supergravity vacuum 
was given in [42,43]. The finite temperature black holes in this asymptotic geometry have 
not yet been found. Since finite temperature breaks supersymmetry and the geometries are 
of cohomogeneity two this is likely to be solved only as a numerical problem. In order to 
ensure supergravity is valid we expect that one will require temperature TX^'^ ^ 1 to avoid 
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a' corrections near the horizon (as in the usual p = MSYM case), and in addition also 
/iA^'^ <^ 1 for the same reason. Thus there will be a one parameter family of solutions 
characterised by the dimensionless ratio x = n/T. 

In the supergravity limit for zero mass fi we know that the energy (= energy density for 
p = 0) and the scalar vev scale as, 



oiV^tWSAVs ^ y^TrW^ ^ bt^/'X'/^ (4.2) 



for t <^ 1 and constants a and b. Adding the relevant BMN mass deformation, the question is 
now how these behaviours will be modified. Suppose we fix the dimensionless ratio x = fi/T 
and vary the temperature. We might expect energy and the scalar vevs again to have some 
power law scaling at low temperature. A priori one would expect that this power law would 
depend on the value of x, so that; 

e = A{x)N^tP'-'^h^/^ , 



\ 7=1,2,3 \ 



E {^^)=B2ix)t<^h'/' (4.3) 



N 

1=4. ..,9 



again for t = T/X^/^, in the limit T , /x < A^/^, so x ~ 0(1), where A{0) = a, Bi{0) = ^2(0) = 
b and p{0) = 14/5, q{0) = r(0) = 2/5. We know of no argument from supergravity that would 
dictate the x dependence of the leading powers p{x), q{x),r{x) in this low temperature limit. 
Now consider the effect of the above mass deformation on our effective action for the 
moduli fields. The cubic Myers term vanishes on the classical moduli action since the $^ are 
diagonal matrices. The bosonic mass terms then generate masses of order ^ for the various 
moduli fields 0^ (again breaking the 5*0(9) global symmetry), which we estimated as. 



^E,classical ^ l j ST^ I ^ a H fi X i ^ Y^ rj.I\2 , f^ ST^ t±I\2 , ^ 



a \ 7=1,2,3 1=4,. ..,9 

ignoring the gauge fields. Now the requirement that x ~ 0(1), i.e. /?/U ~ 1, precisely means 
these two terms are of the same parametric order. Hence the estimate for the classical term 
is unchanged after the addition of the BMN mass. 

We see that strictly the 'moduli' fields (pa are no longer massless moduli, but we shall 
still refer to them as such as we shall shortly see they are parametrically light compared to 
the massive off diagonal modes. Now consider the 1-loop terms. These arise from integrating 
over the heavy off diagonal modes of the matrices. There are now two mass scales, that going 
as \(pa — 4'b\ and given by the moduli, and that arising from the explicit mass /U. The heavier 
of the two mass scales will control the leading behaviour of the 1-loop term. In our estimates 
above, without mass, we found that for t <C 1 then /3|</>a — (l)b\ ^ 1- Now since x ~ 0(1) we 
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have ///3 ~ 1, and hence we see that {(pa — 4'b\ ^ fJ-- Thus our niassless estimate indicates that 
the mass induced by the moduU separation is parametricahy larger than the exphcit BMN 
mass fi in the low temperature regime of interest t <^ 1. Hence adding in this BMN mass 
should make no difference to the leading estimate for the 1-loop behaviour of the moduli 
fields. 

We therefore conclude that for fixed x = fj,/T ~ 0(1) the power law dependence of energy 
(and free energy) on temperature, and also the dependence of the scalar vevs on temperature 
will be as in the massless case. Hence we expect, 



Tt<^i<^i 2 1 

\ ^=1,2,3 \ 



2 1 



J2 (— r^)~i?2(x)t^A3 (4.5) 



N 

7=4. ..,9 



in the limit T , ij. <^ X^'^. Comparing to our naive expectation above, we see that our estimate 
equating tree and 1-loop terms has determined p{x) and q{x) as being constant. Since the 
supergravity black holes have not yet been found we may call this a prediction for supergrav- 
ity.5 

5. Computing the 1-loop corrections 

We now detail the calculation to determine the 1-loop corrections to the classical moduli 
action, the relevant results of which we stated in equations (3.3) and (3.5). This calculation 
is a straightforward generalisation of the calculations considered previously in the context of 
Matrix theory [9,38,46], the principle difference being that our moduli are fields rather than 
being only functions of time, and that we will be interested in calculating both thermal and 
non-thermal corrections. 

Rather than perform a full 1-loop calculation, we will use the shortcut of computing 
the 1-loop correction to a simple classical solution, and then use symmetries to deduce the 
general answer. This will not yield the full 1-loop correction but will give us the information 
we require to make the estimates discussed above. Formally our calculation requires p > 2, 
but since we will keep p general and obtain an answer with analytic dependence on p we may 
deduce the result for any p, and in particular the cases p < 3 we are interested in here. 

We consider off-diagonal fiuctuations 6^, SA^^, 5^ about a classical solution of the moduli 
fields a'ijCpa, so, 

<b = a'iSab + 5i^fe , ^ab = $aSab + ^^ab , ^ab = ^^ ab ■ (5.1) 

Since we are not interested in computing the effective action contribution involving derivatives 
of the gauge moduli, we will take these to be constant and only non-trivial in the Euclidean 



^Currently Costa, Penedones, and Santos [52] are numerically constructing these solutions, and it should 
soon be possible to compare these predictions with their numerical results. 
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time direction so that, 







(5.2) 



for constants a^ which encode the holonomy of the N U{1) gauge fields about the time circle. 
Being non-compact there are no holonomy moduli associated to the spatial directions. Since 
the classical scalar moduli equations of motion are the harmonic conditions d^d^(f)a = 
a linear coordinate dependence gives a solution and allows us to consider gradients in the 
scalars. An important point is that as Euclidean time is periodic, we may only choose linear 
dependence in the spatial directions and not in Euclidean time. We restrict all the moduli (f)a 
to have dependence on two spatial directions, which we take to be x^ and x^ . Then we have. 



ba = qa + X^Va + X^Ua 



(5.3) 



where qa, Va and Ua are constants. Technically this implies that p >2. As stated above, since 
we will keep p general, then in the end we will deduce our answer for general p including the 
cases p < 3. We will now split the world- volume spatial coordinates x* = {x^, x^, . . . , x^} into 



X 



and x^ and the remainder x* = {x^, . . . , x^}. 



Due to the gauge invariance we must introduce Fadeev-Popov ghost fields and a gauge 
fixing term. Following [38] we choose the standard background field gauge choice for the 
fluctuations, 



G = d^'5Af,-i[(l)^,6<!>^] =0 



(5.4) 



adding -^ J dTdx^Tr^G"^ to the action. Having removed the gauge zero modes, the action 
may be expanded in the fluctuations, the leading term being that from the classical moduli 
action so. 



S 



E , classical 



N 



drdx^ ^ y 

a 



(5.5) 



There is no linear term since we are expanding about a classical solution, and at quadratic 
order we find for our bosonic fields, including the gauge fixing term. 



oE 
'-'B,quad 



IN 
2X 



drdxP Y, 



M° 



6a' 



6A^ 



ba^ 



5A^ 



M 



/A„5 

^ab 





V 2i{vi-vi)2i{ui-ui) 






Kb 




6a' 










5A 



ba^ 



5^^ 



ba 



•M- 



Aa6% 





-2i{vi - I 
-2i{ui - I 

Kbhj 



(^Kb\ 

^Kb 
^Kb 
^Kb 



(5.6) 
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with, 

Aa6 = -d^d^ + 2i {aa - at) dr + {aa - a^f + |(?a 

For the fermions we obtain, 

IN 



(5.7) 



^F,quad 



2 A 



JdTdxPY,^^ba{l'^d^ 



i7 " {aa - Ofc) - i7 • Ua - '/'b <J*afe (5.8) 



where the Euchdean gamma matrices ^y'^ are related to the Lorentzian ones as 7'^ = ^7* and 
yi — jt^ g^j^^ finally we have ghost fields, 



'^ghost,quad \ 



/ drd^x ^ CfeaA, 



ib(^ab • 



(5.9) 



Integrating out the off-diagonal fluctuations at quadratic order gives one-loop determinants 
which yield an effective action for the moduli. Now consider for each colour pair {ab) the 
following 2x2 matrix and its eigenvalues Ajafej 



{Va - Vb) ■ {Va - Vb) {Va - Vb) • (Ua - Ub) 

{Va - Vb) ■ {Ua - Ub) {Ua " Ub) ■ {Ua " Ub) 



■ Gi,ab = K,ab Gi,ab , i = 1,2 (5.10) 



Then the bosonic, fermionic and ghost fluctuations result in effective actions which can be 
written in terms of these eigenvalues as, 



oE 
'-'B,quad 



a<b 



i,ab 



'^F,quad — + /^ Tr 
a<b 



6 In (Aab) + J^ In (A^j, + 2y^b) + E ^^ (^«^ " ^A^^ 

j=l,2 j=l,2 

21n(Aa5-F y/Xl^ab + a/A^) -F21n(Aaf,+ v^Ai,a6 - \/A2,a6 



+2 In i^Aab - y^Xi,ab + y%^bj + 2 In (^Aab - \A^ - y^A^ 
^G,,nad = +ETr[21nA 



^afe 



(5.11) 



a<fe 



respectively, and the trace is taken over momenta, with periodic boundary conditions in time 
for the bosons and ghosts, and anti-periodic for the fermions. As detailed in Appendix A 
these traces may be evaluated to yield the 1-loop effective action which consists of colour 
pairwise interactions, 



■^ a<b 



,Va - Vb,Ua - Ub\ 



(5.12) 



15 



where keeping terms involving up to 4 powers in the gradients, Va and Ua, then, 

l-up + \u\ 



L [a, q, V, u] 



r (^^) (2 (\v\'' + |n|4 + 2 (v.u^^^ ^1-12 , u-rl2^2 

., ^ i±p 
(47r) 2 



-/9|<^l 



+ - 



COS (/3a) / /3 



180 



32-^(|.f.|fiP)(A 



S (("■*)' + I" 



+ |a|' + 2(..fi)=) + l(|.„1^ + |„1=f)(^]' 



+ — I I"?! (^•0) +1^1 [u-4>] +2v ■ u [v ■ (f>] [u ■ (p 



P 



•I 2 , \-r\2 
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+ 1^1') 



V- 6] + [u- 






+ ... 



36vv"--^y'+("--^') j g 



=(j+xli/+a;2ji 



where the . . . represent either terms of higher power than four in velocity, or terms that 
are suppressed relative to those displayed in the limit /3|0| » 1 where (/> = q + x^v + x'^u. 
We see that the first term going schematically as ~ {v'^,v'^u'^,u'^)/\(p\'^~'P does not explicitly 
depend on the temperature f3~^ - this is a non-thermal contribution - whereas the remaining 
ones explicitly have temperature dependence - the thermal terms - and are exponentially 
suppressed in the limit fi\4>a — (pbl ^ Ij which from the estimates of section 3 we understand 
as the low temperature limit. 

We have computed the classical and leading one loop corrections up to 4 powers in 
gradients Va, Ua in the low temperature limit f3\(j)a — 4>b\ ^ 1 for a particular classical solution. 
We may use this result to obtain terms in the classical moduli action, and the leading 1-loop 
correction, for a general moduli configuration a{x^), (f){x^) using the fact that, Va = dicpa and 
Ua = d24>a and that spatial rotational SO{p) invariance implies. 



^N 2^ 



/3 



(5.13) 



\v\ [v 



+ \u\ [u 



|^|2 + U|2 = diS-d'c 



(v-$) 



+ [u 



di' 



-\2 



\v\' + \u\' + 2{v.uY = [di(t>-dj 



■ff 



(5.14) 



+ 2v ■ u[v ■ (j)] [u- (f) 



cl)-dict>] U-dj(t>] [ff<p-d^(t> 



Without taking linear dependence on two coordinates (and hence technically requiring that 
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p > 2), but taking only the simpler case of dependence on one, then one would not be able 
to distinguish all the tensor structures above. 

Notice that the derivatives above are spatial di rather than spacetime d^ derivatives. An 
important point is that we did not consider linear r dependence (due to periodic Euclidean 
time), and because the thermal boundary conditions explicitly break the time-space sym- 
metry, we cannot simply deduce the structure of terms involving r derivatives from those 
involving spatial ones, ie. the effective action is not Lorentz invariant. An important excep- 
tion to this is the non-thermal term, which does not explicitly depend on the thermal circle 
period, and arises also for non-compact Euclidean time. In this case the term is not sensitive 
to the thermal boundary conditions breaking time-space symmetry and should exhibit the 
worldvolume Euclidean Lorentz symmetry, and we use this fact to deduce the expression for 
Lnon-thermai bclow in equation (5.17). We also emphasise that the thermal potential term 
which involves no derivatives is fully determined by our analysis. However, generally for the 
thermal terms it is only terms involving spatial derivatives that we can compute, and our 
analysis misses those with time derivatives. 

A second important point is that we only computed terms involving single spatial deriva- 
tives acting on the scalars. Since we computed these terms by considering fluctuations about 
an exact solution with the scalars having linear dependence on the spatial coordinates then 
terms such as didj4>a vanish and hence are not visible to our analysis (obviously unless they 
are related to our terms by integration by parts). We note that at four derivative order the 
non-thermal term is highly constrained by supersymmetry [53-56] and we do not expect such 
terms involving second derivatives. However, we certainly do expect that for the thermal 
contribution there will be terms of such a form that we have missed. 

Hence we may deduce a covariant form for the 1-loop correction to the moduli action 
which fully determines the non-thermal terms but only gives certain contributions to the 
thermal corrections, and ignores derivatives in the gauge fields. We find, 



(5.15) 



where the . . . above indicates we are keeping only the leading terms in the limit /3\(f)a — 4'b\ ^ 1) 
and the Ua are holonomies about the time circle. 



Ua = e'^'^^"- = e^i^'^^'^'^'"' (5.16) 
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and, 



^non— thermal 



(47r) 2 



5„(/>-a^0 a^<^-9> 



du 



^thermal 



PA 



D 



1 



/? 



(27r)2/3 Vl<Al 



|0|7-P 



101 



7-p 



+ higher derivatives 



/? 




^ K- a.*) (8'* ■ 8>$) ( ^) + 1 (a.* . a>f (^) 
-1^ («•*■«•*)(*■»'*)(*■**) (ill) 



1 
36 



• dr 



• • d^(j)\\ -^ + higher derivatives 



(5.17) 



where again we emphasise the derivatives di in the thermal terms given involve only the 
spatial directions, and not time, whereas we expect the non-thermal term to be (Euclidean) 
Lorentz invariant, hence its spacetime derivatives 9^ above. 

It would be interesting to compute the 1-loop thermal corrections about a general moduli 
configuration, rather than the particular one we have chosen, to confirm that the thermal 
terms we have missed, those involving second or higher derivatives, and those involving time 
derivatives, are indeed exponentially suppressed as we assume here. For completeness it would 
be satisfying to also compute the derivative terms in the gauge fields. 

6. Summary 

This paper has been concerned with (1 + j»)-dimensional MSYM at finite temperature in the 
't Hooft limit, where the dimensionless temperature t = T /\^'^^~^> is taken to be finite as 
A^ — 7- oo. For low temperatures t ^ 1 Maldacena duality conjectures the thermal physics is 
described by IIA or IIB supergravity black holes. We have discussed the low energy moduli 
of this theory and their 1-loop quantum corrections, focussing on the scalar moduli. Having 
computed these 1-loop corrections we then argue, following the logic applied for p = in the 
Matrix theory limit [39-41] and more recently by Smilga in the 't Hooft limit we consider 
here [37], that a new scale emerges when this loop expansion becomes strongly coupled, 
which we estimate by equating the classical and 1-loop terms in the moduli effective action. 
Under reasonable assumptions these estimates yield precisely the correct temperature and 
A^ power law behaviour for both the energy density and scalar vevs that are predicted by 
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the dual supergravity black holes. Our canonical approach naturally reproduces both the A^ 
and t dependence, whereas the p = Hamiltonian analysis in [37] assumed the correct A^ 
dependence in deriving the behaviour of t. 

The simple power law predictions of supergravity only apply when it is valid, in the low 
temperature regime where t <^ 1. Our canonical approach allows us to qualify the 1-loop 
corrections into non-thermal and thermal components. We find that the simple scalings pre- 
dicted by supergravity arise from equating the non-thermal corrections (those that arise in the 
zero temperature limit) of lowest derivative order at 1-loop with the classical moduli action. 
Self consistency of ignoring the thermal corrections and higher derivative non-thermal terms 
requires that one is in the low temperature regime, so the dimensionless temperature t <C 1. 
Thus not only are we able to reproduce the supergravity predictions, but we also reproduce 
their regime of validity from the picture of having strongly coupled 1-loop corrections. For 
t ~ 1 the thermal and higher derivative non-thermal corrections would all contribute equally, 
and a much more complicated dependence of energy density and scalar vevs on temperature 
would result, rather than simple power laws. 

The arguments we make are based on simple estimates and are therefore imprecise in 
nature. However, since for all p < 3 the supergravity predictions for the behaviour of both 
energy density and scalars with N and t are reproduced correctly and only in the t <C 1 
regime where they should apply, we are compelled to believe this analysis gives the correct 
account for the origin of the various power law scalings of these physical quantities, and it is 
not simply numerology. 

We have provided an application of these ideas to deduce constraints on the thermody- 
namic behaviour of the IIA supergravity black holes dual to the BMN quantum mechanics 
in the 't Hooft regime. The BMN theory is a mass deformation of the p = MSYM. The 
analysis shows that for fixed ratio of mass to temperature, then the power law scaling of the 
free energy and scalar vevs with varying temperature in the low temperature limit should be 
the same as that of the undeformed p = case (the BFSS theory). Naively one might have 
imagined that this power law would depend on the value of the fixed dimensionless ratio of 
mass to temperature, whereas our analysis indicates that it does not. The dual black holes 
are metrics of cohomogeneity two, and hence will require numerical methods (such as those 
reviewed in [57]) to find them. We understand such a numerical study will soon be complete 
and it will be interesting to see whether our prediction is confirmed [52]. If it is, then an 
obvious question is how this constraint on the thermodynamic behaviour can be understood 
from the supergravity directly. 

Obviously one important open question is to confirm that the scalings we see from 
equating the classical and 1-loop non-thermal terms are compatible with higher loop or- 
ders. Smilga [37] has pointed out that for p = the higher loop corrections [47, 58] appear to 
be compatible with the scaling. We expect that the same will be true for all p. It may be that 
the underlying principle that controls this emergent strong coupling scale is the generalised 
conformal symmetry of Jevicki and Yoneya [59]. Given that for p = 3 it is the conformal 
symmetry that dictates the power law scaling of all quantities with temperature, it seems 
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reasonable that this generahsed conformal symmetry is reahy the structure that controls the 
non-trivial temperature scaling of quantities for p < 3. Thus understanding the role of this 
generalised conformal symmetry seems likely to be a fruitful direction to pursue. 

As outlined in the introduction, the motivation for this work is to provide an avenue to 
improve direct computations for MSYM in order to simulate quantum gravity. Our analysis 
gives insight into the scales in the low energy dynamics relevant for the dual IIA or IIB gravity 
regimes involving black holes, and we are hopeful that this new information - in particular 
the estimates indicating the behaviour of the scalars and their derivatives with temperature - 
may allow improvements in the direct solution methods employed. In the analytic context of 
the mean field theory approach of Kabat, Lifschytz and Lowe [5-8] it may inform improved 
Gaussian approximation schemes. In the Monte Carlo approach [10-13] it may allow resources 
to be focussed more efficiently on the most relevant field configurations. 
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A. Appendix: Evaluating the 1-loop determinants 

At 1-loop for each colour pair we must evaluate a sum of terms involving traces over the 
discrete momentum modes in periodic Euclidean time, and the continuum of spatial momenta, 
that each take the form, 



D{M) = Tt In ( -d'^du + 2iadr + a^ + |(A|^ + M 



(A.l) 



for constants o, g, v, u and M, and, (p = q + x^v + x^u. Then for a boson (periodic in 
Euclidean time) and a fermion (anti-periodic in Euclidean time) we may expand in momenta 
in time as. 



Dboson{M) = ^Tr' In ( -d'di + 

DfermioniM) = ^ Tr' 



27rn 



P 



a] +\6\' + M 



nS; 



In -d'd., + 



27m 



a] +\(f)\^ + M 



(A.2) 



where Tr' represents the remaining trace over the spatial momenta. Analogously to the 
eigensystem defined in (5.10) we define. 



V • V V • u 



V • u u • u 



'^i^i 1 ^ — r , Z 



(A.3) 



where the eigenvectors e^ corresponding to the (real) eigenvalues Ai^2 are chosen to be or- 
thonormal (since the matrix is real symmetric). Using this orthonormal basis we define 
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quantities 01^2, ^1,2 and new coordinates y^'^ by, 



q ■ V 
q- u 



aiei + 0262 



y^ei+y^e2, Ai^2 



T / ^ ■ 



A 



1,2 



n,2 



U- (f> 



.(A.4) 



With these quantities the trace over continuous spatial momenta with the expUcit dependence 
on x^ and x^ may be written using the proper time representation as (see for example [38] 
where there is dependence on one coordinate), 



DbosoniM) = y^ [ cPx I 



4VA1A2 



da / 
(J V (47r)P(TP-2sinh2VArcrsinh2VA^o- 



(A.5) 



X exp 



-a\M+['^-a] + 



^^/tanW5._ 



and for DjermioniM) we have the same with ^^^^ —5- S„ei+z- 

For each colour pair, writing A = —8^8^ + 2iadT + a^ + |g + x-^w + x^-up then the boson 
plus ghost effective action, D '^ , takes the form. 



D 



B+G 



-Tr 



41n(A)+ Y^ ln(A + 2v^j +ln(A-2yA 



1=1,2 



(A.6) 



4i^6oson(0) — DbosonC^y ^1) — ^6oson( — 2a/ Ai) — Dhosoni'^y ^2) — Dhoson{ — '^y ^2 






da 



4VA1A2 



a Y {AT:)PaP~'^sm\).2^/\[asm\i2^/X2a 



4 + 2 cosh 2^X70" + 2 cosh 2i/A^cr 



X exp 



-a I M+( ^-a| + 



V Vf^^^^^^-llA: 



i=l,2 ^ 



A,;0- 



For the fermions we have 

In ("a + \f\i + VA2 ) + In ( A + VAi - ^^2 ) + In ( A - VAi + \/A2 j + In ( A - VAi - \/A2 

Al — \/A2) + 2Dfermion{ — \/M + V A2) 



Z)^ = 2Tr 



^-'-^ fermionyy ^1 ~r y A2J H~ ^ly fermior 

+2Df^r^ion{ — \ ^1 — v A2 



(A.7) 
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We may expand these for small gradients in powers of v and u, which translates into expanding 
the above expression in powers of the eigenvalues Ai^2 (noting that Ai are 0(1) in the limit 
of small V and u) . Working to the order appropriate for four powers in the gradients v and u 
we obtain, 
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B+G 



nei 



Y^ d^x 



da 
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+ ... 



(A. 



where the . . . represent terms involving higher powers of v and u. Now we perform the a 
integration, use the relations. 



Xi + X2 = v-v + u-u, AiAi + A2^2 = (^- 0) +iu- 
{XifAi + {X2fA2 = v-v(v-A +u-u(u-(f\ +2v-u(v-(f))(u-(t) 
(Ai - A2)^ = 2 ({v ■ vf + {u ■ uf + 2 (v.uf) - {v ■ v + u ■ uf 
2X1X2 = — ({v ■ v) + {u ■ u) +2[v.u) \ + {v-v + u-uy 



(A.9) 



and define, D^ = D^^*^, D = D^ , together with. 
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to yield, 
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"90 
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16 

± — \ V ■ V [V 
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+ u ■ u (u ■ (f)] +2v ■ u (v ■ (f)] (u ■ (j)] ] H^ 



^^-{v-v + u-u)i(v-(p] + [u 



Ht 



Tg 1 \V-(t>] +[U 



Ht + ... 



(A.ll) 



As detailed in the Appendix B the infinite sums H^ can be written as, 
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T(q 



An 



l~2q+p 



n=l 



.E_i 



n 



^i-,-f 1^^ 



(A.12) 
cos {nj3a) . 



so that in the limit (3\(j)\ ^ 1 we may expand the modified Bessel function to obtain. 
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and so using f dr = [3 we arrive at, 



j^B+G ^ j^F^j^+^j^ 



dT / d^ 
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36 ' '^" "^ 



/3 



+ 



where we have worked to fourth power in the gradients v, u and at each order in these 
gradients the terms given by the dots . . . are subleading in the hmit that /3|0| ^ 1. 






B. Appendix: Evaluating an infinite sum 

We now consider the infinite sum, 



m = E 



1 



:a (n - &) + M2 



(B.l) 



Since this is periodic in 6 ~ 6 + 1 then we may express it as a Fourier sum, 

+ 00 



F(h)= Y, e^^'^^Fk 



(B.2) 
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with Fourier coefficients, 
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which evaluate to give, 



rl-2«;r(u;-^) 



Fo = V^^ r(^) (B-4) 



and 



-I 



^fc = ^-/^ = ^(^) i^i_^(2vrM|/c|) , A:/0. (B.5) 



Hence, 



F(6) = VvrM--':^^ + 2^^ "' i^._ (2vrM|fc|) cos (2vrfc6) (B.6) 



fc=i 
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